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i-rt ! The evolution of enzyme production is studied analytically using ideas of the 

P5 ' group selection theory for the evolution of altruistic behavior. In particular, we 

O I argue that the mathematical formulation of Wilson's structured deme model ( The 

Evolution of Populations and Communities, Benjamin/ Cumings, Menlo Park, 1980) 
is a mean-field approach in which the actual environment that a particular individ- 
ual experiences is replaced by an average environment. That formalism is further 
J^ I developed so as to avoid the mean-field approximation and then applied to the 

vQ ■ problem of enzyme production in the prebiotic context, where the enzyme pro- 

ducer molecules play the altruists role while the molecules that benefit from the 
0^ ' catalyst without paying its production cost play the non-altruists role. The effects 

^^ I of synergism (i.e., division of labor) as well as of mutations are also considered and 

the results of the equilibrium analysis are summarized in phase diagrams show- 
ing the regions of the space of parameters where the altruistic, non-altruistic and 
j^ ■ the coexistence regimes are stable. In general, those regions are delimitated by 

discontinuous transition lines which end at critical points. 
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1 Introduction 

The controversial issue of the evolution and maintenance of altruism has proba- 
bly entered the field of prebiotic evolution when Maynard Smith [Q remarked that 
giving catalytic support in a molecular catalytic feedback network, such as the hy- 
percycle y], is in fact an altruistic behavior. As a result, such systems are extremely 
vulnerable to the presence of parasites, i.e., molecules that receive catalytic sup- 
port but do not give support to any other molecule in the network. However, the 
stability of this type of cooperative networks is crucial for the theories on the origin 
of life, as Eigen has shown that the lengths of competing self-replicating molecules 
are limited by their replication accuracies and so they cannot integrate sufficient 
information to code for a complex metabolism ^ ^. For the sake of concreteness, 
we define an altruistic behavior as one that is detrimental to the fitness of the indi- 
vidual who expresses it, but that confers an advantage on the group of which that 
individual is a member H. 

In the traditional group selection modeling, based on the Island models of 
Wright P , it is assumed that the population is divided into reproductively isolated 
subpopulations or denies B . The stability of the altruists is achieved by postulating 
the existence of an external extinction mechanism acting on the denies that takes 
place at a rate depending on the deme composition. Of course, such extinctions 
will favor the occurrence of individuals that lower the probability of extinction of 
the deme they belong to which, in the case, are the altruistic individuals R? Hi HI- ^ 
more modern formulation of group selection put forward by Wilson [ pO[ considers 
the denies as trait groups, in which the actual ecological, biochemical or social 
interactions occur, but the individuals are allowed to access and compete for the 
total resources available in the environment. Clearly, in this formulation the notion 
of group or deme is somewhat blurred since, as will become clear in the examples 
given below, there is a stage of the life cycle of the individuals when they leave 
their demes to (effectively) interact with each other. 

Actually, it is not so hard to envision physical systems described by Wilson's 
trait group or structured deme model. For instance, some basic features of viral 
selection dynamics can be modelled by viewing the cells as demes |1^, ^. In this 
case, it is assumed that only N free viruses entry and hence infect a cell; however 
inside the cell the viruses undergo exponential growth leading to the burst of the 
cell and the consequent release of free viruses which will again infect (colonize) 
other cells, and so on. As only N viruses can infect each cell, there is an effective 
competition between all individuals in the population. This restriction, though very 
far-fetched, does not seem to change qualitatively the behavior of the system [pT[ 
and, in addition, it suits very well to describing in vitro serial passage of viruses |12| . 
Other interesting application of Wilson's formalism, which will be the main concern 
of this paper, is the evolution of enzyme production in the prebiotic context p3| , H . 
Here the demes are rock crevices or suspended water droplets of some fixed size. As 
before, although the macromolecules inside the demes undergo exponential growth, 
they are regularly washed away by tides or distributed by winds, and only a small 



fraction of them is then re-adsorbed to the cracks or droplets. Both examples show 
that the spatial localization of viruses or macromolecules facilitates the selection 
against parasites. Henceforth we will refer to the individuals that do not display 
altruistic behavior as non-altruists instead of parasites, since they can subsist even 
in the complete absence of altruists. 

The mathematical formulation of Wilson's structured deme model is centered 
on the concept of the average subjective frequencies of altruists, which are defined 
as the frequencies of altruists experienced by the average altruist and non-altruist 
in the population ||l^ . These quantities differ from the global frequency of altruists 
because the variance of the distribution of the deme compositions is non-zero, i.e., 
the population is not homogeneous. In particular, the stability of the altruists 
is achieved by assuming that the fitness of both altruists and non-altruists are 
proportional to their subjective frequencies. This formulation may be viewed as a 
sort of mean-field approach in the sense that the fitness of a given individual, say a 
non-altruist, in a particular deme is not proportional to the frequency of altruists it 
actually experiences (i.e., the fraction of altruists in its deme) but to the frequency 
of altruists experienced by the average non-altruist in the population. In this paper 
we show that going beyond this mean-field approach does not make the theory 
any more complicated and, in addition, it allows the identification of a recently 



proposed model for the evolution of altruism [g^ [16|, as well as of a population 
genetics formulation of Eigen's model of molecular evolution |1^, as variants of 
Wilson's group selection model. 

The remainder of the paper is organized as follows. In Sec. we present the 
general formalism that takes into account that the fitness of an individual depends 
on the fraction of altruists it actually experiences in its deme. Otherwise the model 
conforms to Wilson's trait group model with the unlimited growth inside the demes 
followed by the destruction of the demes, and the random sampling of N individuals 
to form each new deme. The formalism is then applied to the detailed study of a 
model for the evolution of enzyme production proposed by Michod iQ in Sec. 
|[ Building on the work of Donato |15, Q, in Sec. ^ we apply our formalism to 



investigate the effects of synergism or division of labor in the prebiotic problem 
of enzyme evolution. A variant of the quasispecies model of molecular evolution 
in which the replicating entities are the individual monomers that build up the 
molecules is considered in Sec. o. Finally, some concluding remarks are presented 
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Sec. |. 



2 The model 

The population is composed of an infinite number of demes, each of which is com- 
posed of N haploid, asexually reproducing individuals. The individuals can be of 
two types, A or B, depending on whether they present altruistic or non-altruistic 
behavior, respectively. By definition, altruistic individuals increase the fitness or 
reproductive rate of all individuals in the deme they belong to, but pay a price for 



that by reducing their own fitness. Thus the key ingredient of any group selection 
model is that the fitness of the individuals depend on the composition of the demes, 
which are classified according to the number of altruists they have: there are A^ + 1 
different types, labeled by the integers z = 0, 1, . . . , iV. Hence, an altruistic individ- 
ual living in a deme of type i has fitness FA{i) while a non-altruistic individual living 
in the same deme has fitness Fb(i), with Fsii) > -Fa(«)- Clearly, either in the viral 
dynamics or in the enzyme production problem mentioned before, the occurrence 
of errors in the replication of the individuals (viruses or macromolecules) may have 
important implications to the equilibrium composition of the population. In order 
to take this possibility into account we introduce the mutation rate u € [0, 1/2], 
which gives the probability that type A mutates to type B and vice-versa. 

To derive a recursion equation for the frequency of altruists pt in the population 
at generation t it is more convenient to introduce the frequency of demes with 
i — 0, . . . ,N altruists in generation i, denoted by Yt (i). According to the discussed 
above, and assuming, as usual, non-overlapping generations (i.e., all individuals in 
generation t are replaced by their offspring in generation t + 1) the average number 
of altruists Ma and non-altruists Mb generated during the stage of unlimited growth 
inside the domes are 

N 

Ma = J2 [(1 - ") '^^ (*) +uiN-i)FB (t)] Yt (t) (1) 

i=0 

and 

N 

^B=J2 [(1 - ") (^ - i) Pb (i) + uIFa (i)] Yt (t) (2) 

4=0 

respectively. Hence the global frequency of altruists in the (free) population at 
generation t + 1 is given by 

Ma 



Pt+i = 



where 

N 



Ma+Mb 
1 -2u 



Wt 



Y,iFA{i)Yt{^) (3) 



Wt^Y. [*^^ (*) + (^ - *) Pb (*)] Yt (i) (4) 

i=0 

is the average fitness of the population. The next step in modelling is to distribute 
these individuals in infinite demes, each of which containing exactly N individuals. 
In the absence of additional information, the most conservative assumption that 
can be made about the re-grouping mechanism is that the individuals are picked 
randomly from the (free) population. This leads to the binomial distribution 



N 



\N-i 



yt+ii^ = \ I ) (ft+ir(i-ft+ir"% (5) 



which together with Eqs. (H) and M) allow the complete description of the life cycle 
of the individuals. 

For the sake of completeness and to facilitate comparisons between the two 
formalisms, at this point it is convenient that we introduce the basic ingredients of 
the original structured deme formalism as proposed by Wilson ||l9 . The conditional 
probability distributions of type A given type I ^ A, B at generation t are defined 

by 

which must be interpreted as follows: considering a particular replicator of type I 
then Pt {*|0 is the probability that such a replicator belongs to a deme containing 
i individuals of type A. Hence the average subjective frequency of altruists as seen 
by altruists is given by 

1 ^ 2 

where a^ = X]i*^^t(*) ~ [J2i^^t(i)] is the variance of the deme distribution and 
Pt = J2i'''^tii)/N is the global frequency of altruists in the population. Similarly, 
the average subjective frequency of altruists as seen by non-altruists is 

1 ^ 2 

In the case the demes are assembled randomly, obeying a binomial distribution, 
one has af = Npt{l—pt) so that for large N the subjective frequencies tend 
to the global one. Since fA{t) > Pt ^ fsit), the main point of introducing the 
subjective frequencies is to show that a population structured in groups of distinct 
compositions can simultaneously enhance the effects of the presence of altruists on 
themselves and diminish those beneficial effects on the non-altruists. Of course, the 
assumption that the distribution of deme compositions Yt(j) affects the dynamics 
only through the average subjective frequencies fi{t),l = A, B is too restrictive, 
limiting, for instance, the choices for the dependence of the fitness of the individuals 
on the deme composition. 

3 Evolution of enzyme production 

According to the scenario proposed by Michod iQ , we consider two types of replica- 
tors, A and _B, and assume that only replicator A can produce a catalyst (enzyme) 



which, however, can catalyze the rephcation of both types of repHcators, but with 
different efficiencies. Since rephcator A, which produces the catalyst, must suffer 
some cost in its noncatalyzed self-replication rate, while replicator B attains all the 
benefits of the catalyst without paying the cost for its production, we have here a 
typical situation of altruistic behavior. The cost associated with being altruistic is 
modelled by assigning the noncatalyzed self-replication rate 1 — r, with r G [0, 1], to 
A and the rate 1 to B. Moreover, the rate of catalyzed replication is proportional 
to the concentration of enzymes in the deme, which in turn is proportional to the 
concentration of replicators A in that deme. Hence, assuming that self-replication 
and the replication catalyzed by the enzyme are separate processes, the fitness of 
a replicator I = A, B belonging to a deme of type i can be written as 

Fi{i) ^ 1- air + ki— i ^ ai,ai + I, . . . , N - I + ai, (10) 

where ai = I ii I = A and if Z = B. Here the parameters ki represent the 
beneficial effect of enzyme mediated replication. In particular, ks — implies that 
the enzyme is specific for the replicator which produced it, as in the one-membered 
hypercycle B. However, it seems more plausible to assume that the primordial 
enzymes were some kind of general catalysts which would facilitate the replication 
of a wide spectrum of replicators, so in the following we will assume that kA > ks- 
We note that Michod considers the case fc^ — ks and u = only [Q . The recursion 
equation (j^) thus becomes 

Pt + l^U+{l-2u)— — — -—JJ ^n , (11) 

1 + Pt [kB ~r) + {kA ~ ks) [j^Pt (1 - Pt) + Pt\ 

which is identical to that obtained using Wilson's original (mean-field) formulation 
[[l3[ . In fact, we note that the coefficient of fc^ in the numerator of Eq. (|ll|) can be 
written as Pt/yi(i), so that the rate of increase of altruists in the population due to 
the replication catalyzed by the enzymes is proportional to the average subjective 
frequencies of altruists. 

It is instructive to consider first the case where mutations are not allowed 
{u = 0) since the steady-state equation obtained by setting pt+i = Pt = P* can 
be solved analytically in this case. Explicitly, we find three fixed points: p* = 0, 
p* = 1 and 

r- kA/N 
^ {kA-kB){l-l/N)- ^ ^ 

A physically meaningful fixed point must be in the simplex [0, 1] and satisfy the 
standard stability condition 

dpt+i 



< 1. (13) 

=p* 

We find that p* = is stable for fc^/r < N, while p* — I is stable for fc^/r > 



dpt 

■ kA 
1 + (1 — 1/N) ks/r. Interestingly, for fcs/?' > N there is a region where both 



fixed points are unstable and so the stable one is the intermediate fixed point ( |12[ ) 
which corresponds to a regime of coexistence between altruists and non-altruists. 
These distinct regimes are illustrated in Fig. Ill where we show the steady-state 
frequencies p* for two different values of the initial frequency of altruists. We 
note that in the case u = the analysis is considerably simplified as only the 
ratios ki/r, I = A,B matter for the stability of the fixed points. We identify four 
different phases in the steady-state regime: the pure altruistic phase (A) associated 
to the fixed point p* ~ 1; the pure non-altruistic phase (B) associated to the 
fixed point p* = 0; the coexistence phase (C) associated to the fixed point (|l2|); 
and the phase labeled {A) — (B) where both p* = 1 and p* = are stable. In 
this phase the two kinds of replicators compete such that there is an all-or-none 
selection, though the winner is not determined by the fitness only, but also by 
its initial abundance in the population. In fact, the basins of attraction of the 
two stable fixed points are delimited by the intermediate fixed point (jl^). These 
results are conveniently summarized in a phase diagram in the plane (fc^/^, ^b /r) 
as shown in Fig. g(a). We note that the transitions between phases {B) and (C) as 
well as between phases (C) and {A) are continuous, in the sense that p* increases 
continuously as those transition lines are crossed. It is important to note that even 
in the case of completely nonspecific catalysis Ua = ks the altruistic replicators can 
dominate the entire population provided that the condition k^ > rN is satisfied. 

We turn now to the more general case where the mutation rate u is nonzero. 
The obvious complication in this case is that p — and p = 1 are no longer fixed 
points and so, in principle, the phases identified before cannot be unambiguously 
defined. However, the threshold phenomena observed in the dependence of the 
steady-state frequency of type A replicators on the scaled catalyst specificity kA/r 
(see Fig. |3|) indicates that an unique extension of the definitions of phases (A), [B) 
and {A)—{B) is possible indeed, provided that fc^/r is not larger than some critical 
value. As expected, phase (C) disappears since its defining characteristic, namely, 
< p* < 1, occurs for all parameter settings in the case of nonzero mutation 
rates. The rich interplay between the stable fixed points is illustrated in the phase 
diagrams of Fig. g. The prominent feature of those phase diagrams is the existence 
of critical points at which the two discontinuous transition lines intersect and, as a 
result, above which it is no longer possible to distinguish between phases (A) and 
(B). For fixed u, r and N the critical point coordinates (fc^, k'j^) are determined by 
requiring that the three fixed points of the recursion equation (|I^) collapse into a 
single one. Accordingly, in Fig. ^ we show the critical point coordinates as function 
of the mutation rate u. As expected, for u — we find kg/r = N regardless of 
the value of r. Of particular interest is the mutation rate at which /c^ vanishes, 
henceforth denoted by Ug, as it signalizes the disappearance of all traces of the two 
distinct regimes associated to altruistic and non- altruistic behaviors, leading to the 
phase diagram of Fig. @(d). Interestingly, at this value of the mutation rate we find 
kj^/r — 2N/ {N + 1) independently of r. The dependence of Ug on the altruistic 
cost r for several values of the deme sizes is illustrated in Fig. u. 

The importance of the finitude of the deme sizes N to the stabilization of the 
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Figure 1: Steady-state frequency of type A replicators for u = 0, A'^ = 5, and 
(from left to right) kB/r = 0, 2, 4, 6, 8 and 10. The initial frequencies are 
(a) po = .999 and (b) po = 0.001. The first three lines in part (a) collapse 
into a single line in part (b). 
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Figure 2: Phase diagrams for A^ = 5 and r = 0.1 showing the regions of 
stabihty of the different fixed points for (a) u=0, (b) u=0.005, (c) u=0.01 
and (d) u=0.0158. The intersection point touches the coordinate axis at 

kA/r = 5/3. 




Figure 3: Steady-state frequency of type A replicators in the population for 
u = 0.005, N = 5, r = 0.1, and (from left to right) ks/r = 0, 2, 2.9, 6, 8 
and 10. The initial frequencies are (a) po = 1 and (b) po = 0. 
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Figure 4: Coordinates of the critical point (a) fc^/r and (b) A;^/r as functions 
of the mutation rate u for N = 5 and (from left to right) r = 0.1, 0.2, 0.3, 
0.5, 0.8, and 1.0. At the points where k'^/r = we find k'^/r = 5/3. 
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Figure 5: Mutation rate Ue beyond which the discontinuous transitions dis- 
appear as function of the altruistic cost r for (from bottom to top) N = 5, 



10, 20, 50 and oo. 
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altruists can be appreciated by considering the limit N —> oo, which corresponds 
to a homogeneous population, in the case of absence of mutations u = 0. In fact, in 
this case the fixed point p* = is always stable, while p* — 1 becomes stable only 
for kA > T + kBi which is a very uninteresting situation from the point view of the 
evolution of altruism since in this case the effective fitness of an altruist (1 — ?- + fc^i) 
is larger than the fitness of a non- altruist (1 + ks) belonging to the same deme. 

4 Synergism 

A puzzling problem in evolution is the existence of complex structures that are of 
value to the organism only when fully formed |18|| . It might be possible that enzyme 
production has become a reality due to the combined work of several molecules, 
each being responsible for the synthesis of different pieces of the catalyst. This 
situation of division of labor between the altruists, termed synergism, can result 
in highly non-additive fitness interactions. To model this case, we assume that 
the advantage to the deme is accrued only if the number of altruists reaches some 
minimal value. Explicitly, we will assume that only individuals belonging to denies 
composed of « > im, with im = 0, 1, . . . , iV, altruists have their fitness enhanced: 
for such demes all individuals have their fitness increased by the factor 1/ (1 — c) 
with c G [0, 1]. The dependence of the fitness of types A and B on the composition 
of the deme is summarized by the following equation 

p..N_J l-a;?- if i<im /,,n 

^^'''~ \ {l-air)/{l-c) otherwise ^ ' 

where, as before, r e (0, 1) is the cost for being altruistic and ai = 1 if I = A and 
if I = B. The recursion equation (0) is then written as 

Pt+i=u + {l-2u)- -— ^ ^ _jv _,,^^^ r^. (15) 



{l-c)il-rp,)+cJ:l^,Y,i^){l-j^ 

The formalism based on the average subjective frequencies cannot be applied to 
describe this dynamics because of the highly nonlinear dependence of the fitness Fi 
on the number of altruists in the deme. 

Before we proceed with the analysis of the steady-state solutions of recursion 
equation (|l5|), we must note that the fitness assignment summarized in eq. ( |l4[ ) 
was used by Donato ||l^ in an alternative model for the selection of altruistic 
behavior, which, similarly to Wilson's structured deme model, though not explicitly 
acknowledged by that author, has a stage of the life cycle of the individuals when 
they interact with all other individuals in the population. In fact, this must be 
so because in Donato's model the relative fitness of an individual, which is related 
to the number of offspring it generates, is defined as the ratio between the fitness 
of that individual and the fitness of the whole population pq, n6[. However, that 
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model has two other ingredients that differ from Wilson's: (i) the sizes of the demes 
are not fixed a priori, but there is a maximal deme size that once reached leads the 
deme to split in two smaller ones; and (ii) the offspring of the individuals of one 
deme in one generation form one deme in the next generation. These rules were 
motivated by the analogy with social animals which live in groups not too large and 
whose offspring remain in the group of their parents. An interesting outcome of the 
model is the possibility of stable coexistence between altruists and non-altruists 
within a same group, which is in fact the situation observed in nature since the 
altruistic behavior is usually exhibited only by some individuals in the group. This 
result contrasts with that predicted by the Island group selection models, namely, 
that in the absence of mutations there are either fully altruistic {i = N) or fully 
non-altruistic (i — 0) groups only 0, ||, ||. As we will show in the sequel, using 
the fitness assignment of eq. ( [l4|) this coexistence regime is associated to one of the 
stable steady-state solutions of the recursion equation ( |l5| ) . 

As before, we will consider first the simpler case where u = 0. As expected, 
p = 1 and p = are always fixed points and, depending on the values of the control 
parameters im, c and r, there can be either one or two additional fixed points. On 
the one hand, p* — is always stable for i^ > 1, while for i,„ = 1 it becomes 
unstable in the region c > r. In fact, for fixed r a stable fixed point appears at 
c — r, increasing continuously from to 1 as c increases. This behavior signalizes 
the occurrence of a continuous transition from a regime characterized by fully non- 
altruistic demes only (j>* = 0) to a regime where inhomogeneous demes formed of 
both altruistic and non-altruistic individuals are allowed also (0 < p* < 1). On the 
other hand, p* = 1 is always unstable for im < N, while for im ^ N it becomes 
stable in the region c > r. In this case both fixed points p* — I and p* — are 
stable but, as pointed out before, only one of the two types of individual will take 
over the population. In the other cases 1 < im < N the intermediate, stable fixed 
point < p* < 1 appears in a discontinuous manner, i.e., p* is non-zero already at 
the outset. 

In Fig. |6| we present the transition lines separating the region in the plane (c, r) 
where the altruistic individuals persist in the population (region below the curves) 
from the region where the only stable fixed point is the non-altruistic one p* = 0. 
Since those curves satisfy c > r, it seems that the surviving altruists are those 
belonging to demes with i > im since they have a larger fitness than non-altruists 
living in demes with i < im- Interestingly, the size of the region of existence of 
altruists decreases with increasing im, reaches a minimal value for im — N/2, and 
then increases towards its initial size as im approaches N (the transition lines for 
im = ^ and im = N coincide). However, it must be noted that the most favorable 
situation to the altruists is the case of no synergism z^ = 1, since only then the fixed 
point p* = 0, associated to the non-altruistic regime, becomes unstable. Moreover, 
the basin of attraction of the intermediate fixed point decreases with increasing 
im and so, unless there is already a large number of altruists at the outset, the 
non-altruists will take over the population. For instance, for «,„ = N the basin of 
attraction of p* = 1 is vanishingly small close to the transition line c — r. This 
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Figure 6: Transition lines for N = 20, u = and, from top to bottom, 
im = 1, 2, 4, 10 (solid lines), and im = 19, 17 (dashed lines). The curves 
for ijn = 20 and im = ^^ coincide with those for ?„ = 1 and im = 10, 
respectively. 
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Figure 7: Phase-diagrams for A^ = 20 and im 
stability of the different fixed points for (a) u 
u = 0.0282 and (d) u = 0.0330. 
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rather frustrating result simply reflects the difficulty already pointed out in the 
beginning of the section, of evolving a synergistic system in nature. A possible 
solution to this problem is provided by the so-called Baldwin effect [|l^ which, in 
the framework proposed by Hinton & Nowlan M, 



20 , assumes the existence of 



a third type of individual, say X, which either by learning, guessing or imitation 
can act as an individual of type A ox B but whose offspring are, of course, of type 
X. These plastic individuals may provide the appropriate conditions (i.e., a large 
number of altruistically behaving individuals) to start the synergistic effects and, 
once this is done, they will become extinct due to the competition with the born 
altruists, leaving thus no trace of their early presence in the population. We will 
leave the investigation of this avenue of research for a future contribution. 

Taking into account the effect of mutation [u > 0) leads to a very rich in- 
terplay between the different steady-state regimes of the recursion equation (nsh 
as illustrated by the phase-diagrams shown in Fig. 1^. In the absence of mutation 
the phase labeled {B) is associated to the non-altruistic regime characterized by 
the fixed point p* = 0, while phase (C) is associated to the coexistence regime 
characterized by the intermediate fixed point < p* < 1. As before, although 
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Figure 8: Mutation rate Ue beyond which the discontinuous transitions dis- 
appear as function of im ioi N = 5 (A), N = 10 (y), N = 20 (□), and 
N = 30(0). The hues are guides to the eye. 



for non-zero mutation rates p = is no longer a fixed point it is still possible to 
distinguish between the fixed points corresponding to the non-altruistic and the 
coexistence regimes, due to the occurrence of threshold phenomena similar to those 
shown in Fig. 0. The main effect of mutation is to produce, at the expense of 
phase (C), a bounded region, labeled {B)~{C), where both phases are stable. This 
region is delimited by two discontinuous transition lines that intersect and end at 
two critical points. As the mutation rate u increases the size of the bounded region 
decreases and disappears altogether at the critical end point Wg at which the two 
critical points coalesce. Hence for m > Me it is no longer possible to distinguish 
between phases {B) and (C). The dependence of Ue on im is shown in Fig. 0. As 
expected Ue = for i™ = 1, regardless of the value of the deme size N , since the 
transition between those two phases is continuous already for u = 0. 
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5 Quasispecies model 

Other interesting application of the forniahsm presented in Sec. is the study of 
the error threshold transition in Eigen's molecular quasispecies model |3|] . Such ap- 
proach has recently been proposed as an (uncontrolled) approximation to the origi- 
nal kinetics formulation of the quasispecies model |17| , without the realization of its 
close connection with Wilson's trait group framework. In this case, the monomers 
play the role of the individuals, and the molecules the role of the demes. However, 
there is no distinction between altruistic and non-altruistic monomers (i.e., there 
is no altruistic cost) but the self-replication rates of the monomers depend on the 
molecule they belong to. Thus, contrasting to Eigen's original proposal, in this for- 
mulation the molecules are not self-replicating entities, being only passive carriers 
of monomers. In this context, it is more appropriate to think of the mutation rate 
u as the replication error rate per monomer. Explicitly, for the single-sharp-peak 
replication scenario we have Fb (i) = 1 — s, i = 0, . . . , A^ — 1 and 

,.. / 1 - s if i = 1,.. .,iV- 1 , , 

FA{^) = ^^ if,^^ (16) 

where < s < 1 is the selective advantage of the so-called master molecule, namely, 
the molecule composed of A'' monomers of type A ||j] . The general recursion equa- 
tion (^) then becomes 

...r="+(l-2.)il-il^i±#. (17) 

1 — s -f sp" 

The only stable fixed point for u = is p* = 1 which corresponds to the domination 
of the population by the master molecules. As the mutation rate increases, two 
distinct regimes are observed in the composition of the population: the quasispecies 
phase (Q) characterized by the master molecule and its close (in the sense of the 
Hamming distance) neighbors, and the uniform phase [U) where the 2^ molecules 
appear in the same proportion. More pointedly, phase [Q) is associated to the fixed 
point p* « 1 and phase {U) to p* ~ 1/2. In Fig. g we illustrate the dependence 
of the steady-state molecule frequencies F(j), given in Eq. (||), on the error rate 
u. Although these results show a remarkable similarity to those obtained with 
the original kinetics formulation of the quasispecies model B, the agreement is 
qualitative only: a full analysis of the location of the error-threshold transition for 
the original model indicates that the predictions based on the recursion equation 
( [l7[ ) are very inaccurate pl[ |. 

In the present study the error-threshold transition corresponds to the discon- 
tinuous transition between the phases {Q) — {U) and [U) (see Fig. M). As in the 
previous models, the discontinuous transition lines intersect and end at a critical 
point (u'^, s'^) given by Q 

1 fN + 1"- ^ 



'-2 F^l <'«' 
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Figure 9: Steady-state frequencies of molecules composed of 10 (master), 9, 
8, 7, 6 and 5 monomers of type A as functions of the error rate for iV = 10 
and s = 0.9. The initial type A monomer frequency is pq = 1. 
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Figure 10: Phase diagram for A^ = 10 showing the regions of stability of the 
quasispecies (Q) and uniform (JJ) regimes. The discontinuous transitions 
end at the critical point u'^ = 0.251 and s'^ = 0.983. 



and 



= 1 



■.N 



N -1 

N + 1 



^ (iV - 1)^ - AN 



(19) 



It must be noted, however, that the rich interplay between phases (Q) and (U) 
depicted in the phase diagram of Fig. |l^ is a consequence of Wilson's trait group 
framework, in which the molecules are disassembled and then randomly assembled 
during the life cycle of the population. Clearly, the use of that framework is inade- 
quate in the context of the quasispecies model, in which the molecules and not the 
monomers are the self-replicating entities. 

6 Conclusion 

In this paper we have basically attempted to re-interpret and unify several models 
dealing with the evolution of altruistic behavior 0, 11^, ll^ in a single framework, 
namely, the 'extended' Wilson's structured deme model of group selection pO| . 
In doing so, we have carried out a thorough analysis of the steady-state regime 
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of a model for the evolution of enzyme production proposed originally by Michod 
[ p3[ , without resorting to the mean- field approximation implicit in Wilson's concept 
of average subjective frequencies [ pO| . Furthermore, the effect of synergism (i.e., 
division of labor) was considered by assuming that the presence of altruists accrues 
benefits only to groups containing some minimal number of that type of individual, 
following thus Donato's alternative group selection model ||l^, 111]. In particular, we 
have obtained the phase-diagrams showing the regions of stability of the altruistic 
and non-altruistic regimes. A particularly relevant result is the finding of a regime of 
stable coexistence within a same group of altruists and non-altruists which, though 
expected from observation, is not predicted by the Island group selection models 
[0, H, ||]. We have also identified a recently proposed variant of the quasispecies 
model, in which the macromolecules are viewed as vehicles for the self-replicating 



monomers 17 , as a particular realization of the 'extended' Wilson's structured 
deme formalism presented in this paper. In addition, we have found that taking 
into account the possibility of mutations leads to interesting qualitative changes on 
the steady-state regime of the model dynamics as, for instance, the appearance of 
critical points in the phase-diagrams of the models. In particular, we have shown 
that there is a value of the mutation rate Ug (see Figs, ra and ^) above which the 
selective pressures are no longer operative, in the sense that it is no longer possible 
to distinguish between the altruistic and the non-altruistic regimes. 

To conclude, we note that in the prebiotic context error-prone replication (mu- 
tation) has played a crucial role in revealing the limitations of non-cooperative 
molecular systems, such as Eigen's quasispecies model, to function as efficient in- 
formation integrators P, H. Furthermore, it was shown recently that mutation can 
have disastrous effects over the stability of altruistic demes in the more traditional 
Island formulation of group selection theory B. In view of this, mutation should 
not be viewed as merely another complication to be added to a model, but as a 
basic test for probing the robustness of any model of integration of information in 
prebiology, this being thus the main motivation for the present contribution. 
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